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determines  the  eigenvalue  of  each  of  the  selected  component  modes  of  the 
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mode  conversion  approach.  The  resulting  algorithm  provides  a  reliable 
and  cost-effective  prediction  of  amplitude,  phase  and  MIPD  of  an  Omega 
station  signal  for  a  given  location  and  time. 
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1  . 


INTRODUCTION 


The  worldwide  navigation  position- fixing  capability  of 
the  Omega  Navigation  System  is  predicated  on  the  assumption 
that  the  Omega  signal,  a  multimode  signal  (in  the  waveguide¬ 
mode  formulation),  can  be  adequately  described  by  the  Mode  1 
component*  of  the  signal  for  any  location  and  time.  The  lack 
of  Mode  I  dominance  in  a  signal  is  known  as  "modal  interfer¬ 
ence".  As  a  consequence  of  modal  interference,  the  phase  of  a 
received  Omega  signal  can  deviate  (differ)  significantly  from 
the  modeled  phase  (which  is  assumed  to  be  the  phase  of  the  Mode 
1  component)  typically  used  in  computing  an  Omega  position  fix. 
Thus,  in  a  fix  computation  any  use  of  signal (s)  having  exces¬ 
sive'1"  modal  interference  -  induced  phase  deviation  (MIPD)  is  ex¬ 
pected  to  degrade  the  fix  accuracy. 

The  objective  of  the  effort  reported  herein  is  to  de¬ 
velop  a  cost-effective  algorithm  for  predicting  amplitude, 
phase,  and  MIPD  of  an  Omega  station  signal  for  a  given  loca¬ 
tion  and  time.  This  algorithm  is  required  for  generating  a 
worldwide  database  of  the  amplitude  and  MIPD  of  10.2  and  13.6 
kHz  Omega  system  station  signals.  The  database  is  currently 
being  developed  by  The  Analytic  Sciences  Corporation  (TASC) 
for  the  Omega  Navigation  System  Center  (ONSC),  the  U.S.  Coast 


*  Mode  1  is  the  lowest  phase - ve loc ity  transverse -magnet ic 
mode  of  the  multimode  signal  propagating  in  the  "earth- 
ionosphere"  waveguide  formed  along  the  signal  path  (Ref.lj. 

+  A  signal  is  considered  to  have  excessive  MIPD  if  it  exceeds 
the  ONSC - recommended  threshold  value  of  20  centicycles 
(Ref.  2 ) . 
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Guard  organi zat ional  unit  which  is  responsible  for  operating, 
maintaining,  and  improving  the  Omega  system. 

An  overview  of  the  Omega  signal  propagation  mechanism 
including  modal  interference  phenomena  is  given  in  Chapter  2. 
The  signal  prediction  algorithm  including  the  assumptions  and 
approximat ions  used  to  develop  the  algorithm  are  presented  in 
Chapter  3.  The  algorithm  development  effort  is  summarized  in 
Chapter  4 . 
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OMEGA  SIGN'LL  PROPAGATION 


Section  2.1  presents  an  overview  of  the  Omega  signal 
propigation  mechanism.  (A  detailed  discussion  of  the  VLF / 
Omega  propagation  theory /mechanism  can  be  found  in  Refs.  1, 
and  3  through  8.)  A  discussion  of  modal  interference  phenom¬ 
ena  and  the  resulting  signal  propagation  characteristics  is 
provided  in  Section  2.2. 


2.1  SIGNAL  PROPAGATION  MECHANISM 

Omega  signals  propagate  in  the  space  between  the 
earth's  surface  and  the  D-region  of  the  ionosphere,  known  as 
the  "earth- ionosphere"  waveguide  (see  Fig.  2.1-1).  Signal 
propagation  in  the  waveguide  formed  along  the  signal  path  is 
conveniently  described  using  a  full-wave  waveguide-mode  repre¬ 
sentation  for  the  signal.  t n  this  representation,  a  signal  is 
expressed  by  a  sum  of  the  "modes"  of  the  waveguide.  Each 
waveguide  mode,  a  solution  of  the  waveguide-mode  equation,  is 
uniquely  characterized  by  its  complex  eigenvalue.  The  lo¬ 
cally-varying  electromagnetic  properties*  of  the  signal  path 
determine  the  propagation  characteristics  (i.e.,  amplitude  and 
phase  variations)  of  the  signal's  component  modes  and  the  re¬ 
sulting  multimode  signal  along  the  path.  Signal  propagation 
characteristics  of  a  mode  are  determined  by  four  locally-vary¬ 
ing  modal  parameters:  attenuation  rate,  phase  velocity,  exci¬ 
tation  factor,  and  height-gain  function.  Modal  parameters 


-Ground  conductivity,  solar  illumination  condition,  geomag¬ 
netic-field  magnitude  and  orientation  relative  to  the  signal 
propagation  direction. 
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Figure  2.1-1  Path-Waveguide  Geometry 

depend  on  the  eigenvalue  and  the  electromagnetic  properties  o 
the  path' 

Since  an  Omega  signal  path  is  inhomogeneous  due  to 
local iy- varying  properties  of  the  pcth ,  available  homogeneous 
waveguide  solution  techniques  must  be  modified  to  include  the 
effects  of  signal  scattering  occurring  at  the  path  property 
d  iscon.t  inu  it  ies  There  are  two  commonly -used  approaches  for 
computing  the  scattering  effects.  one  uses  a  mode  conversion 


Henceforth  in  this  report  the  words  path,  and  waveguide  a  re 
1 nte rchangeab 1 y  used  to  denote  the  waveguide  formed  along 
t  h  e  s  i  g  n  a  1  p  a  t  h  . 
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approach  (Ref.  9);  while  the  other  employs  a  '.venze  1  -Kramers  - 
Brillouin  (WKB)-type  approximation  (Refs.  1,  8,  and  10  j. 

In  che  mode  conversion  approach,  modes  are  assumed  to 
be  uniquely  and  independently  identifiable  in  a  local  sense. 
Even  though  the  characteristics  of  a  mode  depend  only  on  local 
properties  of  the  path,  the  propagation  of  each  mode  from  one 
point  to  the  next  point  along  a  path  depends  on  the  local  re¬ 
distribution  of  the  signal  energy  occurring  among  its  compo¬ 
nent  modes.  The  redistribution  is  caused  by  the  mode  conver¬ 
sion  (see  Section  2.3)  which  occurs  at  a  path  point  with  rap¬ 
idly-varying  mode  - e lgenvalue ( s ) . 

In  the  WKB - approximat ion ,  a  mode  is  assumed  to  be 
uniquely  and  independently  identifiable  everywhere  along  a 
path.  Furthermore,  each  mode  is  assumed  to  depend  only  on  lo¬ 
cal  properties  of  the  path  and  to  propagate  independently  of 
the  existence  of  any  other  modes.  This  approach  is  applicable 
to  an  Omega  signal  path  with  gradual ly -vary ing  mode-eigen¬ 
values.  Because  scattering  effects  in  the  multimode  signal 
have  to  be  computed  at  each  point  along  a  path,  a  mode  conver¬ 
sion-based  multimode  signal  propagation  prediction  approach  is 
generally  much  more  expensive  than  a  WKB - approx imat i on-based 
approach . 

A  cost-effective  alternative  for  determining  signal 
behavior  along  a  path  is  to  use  a  combination  of  these  two  ap¬ 
proaches  as  described  in  Chapter  3.  The  next  section  presents 
a  review  of  modal  i nto r f e rence  phenomena. 
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2.2  MODAL  INTERFERENCE  PHENOMENA 

Although  an  Omega  station  signal  is,  in  general,  a 
multimode  signal  at  most  locations  outside  the  " near - f ie Id* " 
region  of  the  station,  the  signal  can  be  adequately  approxi¬ 
mated  by  its  Mocie  1  component .  Thus  the  presence  of  higher- 
order  modes  in  the  signal  can  be  effectively  ignored.  Mode  I 
is  the  lowest  phase- velocity ,  transverse-magnetic  mode  of  the 
signal  along  the  path,  and  it  usually  has  the  lowest  signal 
attenuation  rate.  The  phase  of  Mode  1  signals,  and  hence  the 
phase  of  Mode  1 -dominated  signals,  vary  almost  linearly  with 
distance  from  the  station. 

In  the  near-field  region  of  a  station,  as  well  as 
along  nighttime  and  transition  paths  at  certain  azimuths  (es¬ 
pecially  those  paths  emanating  from  the  low-geomagnetic- lati - 
tude  Omega  stations),  Mode  I  is  often  not  the  dominant  signal 
mode  (Ref.  10).  This  lack  of  Mode  1  dominance  is  called  "mo¬ 
dal  interference",  and  the  resulting  "modal ly - d i sturbed "  sig¬ 
nal  is  composed  of  several  competing  modes  which  may  alter¬ 
nately  dominate  on  different  segments  of  the  path.  In  this 
case,  due  to  differing  phase  velocities  of  the  signal's  compo¬ 
nent  modes,  amp  1 i tude /phase  of  the  total  (multimode)  signal 
exhibits  an  oscillatory  behavior  with  distance  (see  mode-sum 
signal  variations  in  Fig.  2.2-1).  Alternatively,  the  modally- 
disturbed  signal  may  be  dominated  by  a  single,  higher-order 
mode  (e.g. ,  Mode  2  in  the  signal  behavior  shown  in  Fig.  2.2-2; 
with  a  1  inear ly- vary ing  phase-versus-distance  relationship 
which  is  usually  significantly  different  from  that  of  Mode  1. 


'  A  station  near- field  region  typically  extends  from  the 
station  onward  to  distances  of  500-1000  km  along  "all-day 
paths,  and  1000-2000  km  along  "all-night"  paths. 
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An  Example  of  Modal  Interference  in  the  Signa 
(Mode  2  is  Dominant  Mode  of  Signal  Beyond  5 
Megameters  from  Transmitter 
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Modal  interference  in  a  signal  propagating  along  a 
path  from  a  transmitt ing  station  to  a  given  point  is  character¬ 
ized  as  spatial  or  temporal,  depending  upon  the  solar  illumina¬ 
tion  condition  along  the  path  between  the  station  and  the 
point.  The  interference  along  a  path  is  classified  as  spatial 
interference  if  the  path  has  an  "all-day"  or  "all-night"  illu¬ 
mination  condition  as  the  magnitude  of  interference  at  each 
path  point  is  nearly  constant  in  time  (and  thus  depends  omy  on 
spatial  coordinates).  If  the  path  is  in  transition  (i.e.,  the 
path  has  a  day/night  terminator  crossing),  the  interference  is 
referred  to  as  temporal  interference.  The  magnitude  of  the 
temporal  interference  at  a  point  will  generally  vary  in  re¬ 
sponse  to  the  movement  of  the  day/night  terminator  along  the 
path  to  the  point. 

Signal  phase  deviations  arising  from  spatial  interfer¬ 
ence  are  generally  larger  in  magnitude  and  persist  to  longer 
distances  along  a  signal  path  during  night  than  during  day. 
Temporal  interference  arises  from  the  mode  conversion  effects 
occurring  at  the  day/night  terminator  discontinuity  along  a 
transition  path  (see  Fig.  2.2-3)  which  is  subject  to  the 
"usual"  spatial  interference  when  fully  dark..  Due  to  mode  con¬ 
version,  the  energy  of  each  incident  mode  at  the  path- termina¬ 
tor  crossing  is  converted  (or  distributed)  among  several 
transmitted  modes  (reflected  modes  are  small  at  Omega  frequen¬ 
cies  and  are  neglected)  which  propagate  beyond  the  terminator. 
Thus,  because  of  mode  conversion,  an  Omega  station  signal  with 
a  single  dominant  mode  propagating  toward  a  terminator  can  ex¬ 
hibit  significant  modal  interference  (i.e.,  presence  of  several 
competing  modes)  in  the  vicinity  of,  and  beyond  the  terminator. 
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Figure  2.2-3  Earth  - 1 onosphere  Waveguide  with  Day/Night 

Terminator 
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3  .  MULTIMODE  SIGNAL  PREDICTION  ALGORITHM 


3.1  INTRODUCTION 


The  electric/magnetic  field  of  an  Omega  signal  is  a 
vector  quantity  composed  of  three  orthogonal  field  components 
Each  component  is  described  by  a  complex  quantity  with  an  as¬ 
sociated  magnitude  (amplitude)  and  phase.  The  algorithm  pre¬ 
sented  herein  is  used  to  generate  a  worldwide  database  of  the 
amplitude  and  MIPD  (modal  interference  -  induced  phase  devia¬ 
tion)  of  the  Omega  system  station  signals.  The  signal  field 
(also  referred  herein  as  simply  the  signal)  of  interest  for 
generating  the  database  is  the  vertical  component  of  the  elec 
trie  field  of  the  signal  at  the  earth's  surface.  Therefore, 
subsequent  algorithm  development  is  focused  on  the  vertical 
component  of  the  electric  field  of  the  signal. 


Since  Omega  signal  paths  are  inhomogeneous,  homogene¬ 
ous  waveguide  models  cannot  be  used  to  determine  signal  behav 
ior  along  a  path.  A  commonly-used  approximation  is  to  parti¬ 
tion  the  inhomogeneous  path  into  a  large  number  of  segments 
each  with  constant  (electromagnetic)  properties  associated 
with  the  mid-point  of  the  segment.  As  a  consequence  of  this 
partitioning  scheme,  path  properties  are  homogeneous  within 
each  segment  but  change  discont inuous ly  at  the  interface  join¬ 
ing  the  segments.  The  multimode  signal  propagation  character 
istics  along  the  partitioned  inhomogeneous  path  is  determined 
using  the  following  three-step  approach.  First,  multimode 
signal  behavior  is  obtained  along  each  of  the  homogeneous  seg¬ 
ments  of  the  path  using  a  homogeneous  waveguide  model.  Sec¬ 
ond,  the  scattering  effects  in  the  signal  occurring  at  the 
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segment  interfaces  along  the  path  are  determined  using  a  mode 
conversion  model  or  a  WKB - approximat ion .  Third,  and  final, 
the  multimode  signal  propagation  characteristics  along  each 
segment  are  modified  to  incorporate  the  scattering  effects  at 
the  segment  interfaces. 

A  conceptual  illustration  of  the  scattering  effects 

due  to  the  mode  conversion  occurrence  at  an  interface  between 

the  neighboring  dissimilar  property  path  segments  (1)  and  (2), 

is  shown  in  Fig.  3.1-1.  The  incident  and  transmitted  signals 

are  assumed  to  be  composed  of  Modes  1  and  2,  only.  In  this 

illustration,  E^  and  E^  (k  =  1,2)  are  the  complex  values  of 

the  incident  (i)  and  transmitted  (t)  mode  signals  associated 

with  kth  mode  at  the  interface;  S-  .  is  the  complex  scatter- 

J  *  f  h 

ing  coefficient  which  is  the  ratio  of  the  j  mode  of  the 

f  S 

transmitted  signal  and  k  mode  of  the  incident  signal. 
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Figure  3.1-1  An  Illustration  of  Mode  Conversion  Phenomenon 

Occurrence  in  a  Two-Mode  Signal  (Reflected 
Modes  Are  Not  Shown) 
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This  three-step  approach  for  determining  multimode 
signal  propagation  characteristics  along  an  inhomogeneous  sig¬ 
nal  path  is  currently  implemented  in  the  "Fast  Mode  Conver¬ 
sion"  ( FASTMC ;  Ref.  9)  model  developed  by  the  Naval  Ocean 
Systems  Center  (NOSC) .  Although  FASTMC  provides  reliable  sig¬ 
nal  predictions,  it  is  computationally  very  expensive. 

A  cc st-effective  alternative  for  computing  scattering 
effects  at  a  path-segment  interface  with  no  noticeable  loss  in 
the  accuracy  of  the  resulting  signal  is  to  proceed  as  follows: 


•  Identify  the  expected  types  of  segment  in¬ 
terfaces  (along  worldwide  Omega  signal 
paths)  where  the  off-diagonal  elements,  S. 

(j  ^  k)  ,  of  the  scattering  coefficient  J ' 
matrix  are  likely  to  be  significant 

•  Compute  scattering  effects  at  a  path  inter¬ 
face  using  a  mode  conversion  model,  similar 
to  the  one  implemented  in  FASTMC  (Ref.  9), 
if  the  scattering  coefficient  elements,  S. 

(j  *  k)  ,  at  the  interface  are  expected  J ' 
to  be  significant;  otherwise,  use  a  WKB  ap¬ 
proximation  similar  to  that  implemented  in 
the  Integrated  Propagation  Prediction  (IPP) 
model  developed  by  NOSC  (Ref.  11). 


The  Multimode  Signal  Prediction  (MSP)  algorithm  devel¬ 
oped  herein  is  basically  the  FASTMC  model  with  the  provision  of 
computing  scattering  effects  using  the  WKB- approximat ion  of  IPP 
(instead  of  FASTMC * s  mode  conversion  algorithm)  whenever,  based 
on  a  priori  information  about  the  scattering  potential  of  the 
interface,  the  off-diagonal  scattering  coefficient  matrix  ele¬ 
ments  are  expected  to  be  negligible. 

This  alternative  procedure  yields  signal  predictions 
within  one  dB  of  the  FASTMC  predictions  at  about  one-tenth  of 
the  FASTMC  executes  time/cost  when  comparison  is  made  over  a 
representative  set  of  worldwide  Omega  station  signal  paths. 
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The  algorithm  is  designed  so  that  repeated  signal  calculations 
along  the  same  path  (e.g.,  at  different  times  of  day/year)  can 
be  executed  for  substantially  less  cost/time  than  comparable 
executions  of  FASTMC.  This  algorithm  will  be  used  to  generate 
a  worldwide  database  of  10.2  and  13.6  kHz  signal  path  calcula¬ 
tions  for  each  of  the  eight  Omega  stations  for  96  times  of  the 
year;  the  total  computer  cost  is  projected  to  be  about  one 
percent  of  the  comparable  FASTMC  cost,  which  represents  a  pro¬ 
jected  savings  of  several  hundred  -  thousands  of  dollars  to 
ONSC  . 


Section  3.2  presents  the  assumptions  and  approxima¬ 
tions  used  in  development  of  the  MSP  algorithm.  The  list  of 
the  path-segment  interface  types  where  scattering  effects  must 
be  computed  using  the  mode  conversion  model  is  given  in  Sec¬ 
tion  3.3.  The  MSP  algorithm  development  is  presented  in  Sec¬ 
tion  3.4. 


3.2  APPROXIMATIONS  AND  ASSUMPTIONS 

The  MSP  algorithm  is  developed  using  the  following 
approximations  and  assumptions  concerning  the  signal  path 
e lectromagnet ic  properties,  and  signal  mode  structure: 


1)  A  day  (night)  path  segment  is  defined  by 

0°<  x*  <  (90°  <  x  ^  180°)  for  all  x  along  the 

path  segment  and  a  transition  segment  is  one  in 
which  x  -  90°  at  some  point  in  the  segment. 

2)  The  day/night  ionosphere  boundary  of  the 
"earth- ionosphere "  waveguide  (formed  along  the 
signal  path)  is  a  vert ical ly - strat i f ied ,  hori¬ 
zontal  ly -homogeneous  (in  the  direction  of  sig¬ 
nal  propagation)  medium  consisting  of  free 
electrons  (in  addition  to  positive  and  negative 
ions)  whose  day/night  density  and  associated 


*  X  denotes  solar  zenith  angle, 
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collision  frequency  vary  exponentially  with 
vertical  height  (above  the  earth's  surface). 
Such  a  day/night  ionosphere  boundary,  in  terms 
of  Wait  and  Spies's  notations  (Ref.  11),  is 
characterized  by  a  day/night  pair  of  parame¬ 
ters:  ionospheric  reflection  height  and  iono¬ 
spheric  conductivity  gradient.  The  day/night 
values  of  these  parameters  used  in  the  algo¬ 
rithm  are  listed  in  Table  3.2-1. 

3)  The  transition,  i.e.,  change  from  day  to  night 
solar  illumination  condition,  along  a  path  is 
assumed  to  be  a  step  change  in  the  ionospheric 
parameters  from  the  day  to  night  values. 

4)  Signal  at  a  day  path  segment  is  comprised  of 
the  first-three  lowest  phase- velocity  modes, 
Modes  1  through  3  (Refs.  13  and  14)  of  the 
earth- ionosphere  waveguide  characterizing  the 
path  segment . 

5)  Signal  at  a  night  path  segment  is  comprised  of 
the  first-five  lowest  phase- velocity  modes. 
Modes  1  through  4  and  Mode  X  (Refs.  12  and  13), 
of  the  earth- ionosphere  wave-guide  characteriz¬ 
ing  the  path  segment. 

6)  Mode  X ' 3  contribution  to  the  multimode  signal 
is  negligible  along  all  path  segments  except 
for  the  nighttime,  westerly-directed  (i.e., 
geomagnetic  bearing*  between  190°  and  350°) 
segments  in  the  equatorial  (i.e.,  between  10°N 
and  10°S  geomagnetic  latitudes)  region,  as  dis¬ 
cussed  below. 


TABLE  3.2-1 

DAY  AND  NIGHT  IONOSPHERIC  PARAMETER  VALUES  USED  IN  THE  MSP  ALGORITHM 


IONOSPHERE 

ILLUMINATION* 

IONOSPHERIC  PARAMETERS  j 

IONOSPHERIC  REFLECTION  HEIGHT 
(km) 

IONOSPHERIC  CONDUCTIVITY  GRADIENT 
( km  ) 

Day 

(0°<  x  <  90° ) 

70 

00 

o 

Night 

1 90°  <  x  <  180°) 

87 

0.5 

Trans i t ion 

1 X  =  90°) 

Step  change  from  70  to  87 

Step  change  from  0.3  to  0.5 

*  x  denotes  solar  zenith  angle. 


*  Geomagnetic  bearing  is  measured  clockwise  from  the  geomag¬ 
netic  north. 
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In  the  waveguide-mode  theory  representation,  the  odd-  and 
even-numbered  modes  correspond  to  the  transverse -magnet ic  (TM) 
and  transverse-electric  (TE)  modes,  respectively,  of  the 
earth- ionosphere  waveguide.  The  total  signal  along  a  path 
segment  is  commonly  expressed  in  terms  of  the  TM  and  TE  modes 
of  the  waveguide  modeling  the  path  segment.  Mode  X  is  a  new 
mode  recently  identified  by  TASC  and  it  is  found  to  exist  only 
along  a  nighttime  path  and  not  along  a  daytime  path  (Refs.  12 
and  13).  The  dependence  of  the  Mode  X  eigenvalue  and  modal 
parameters  on  the  path  properties  is  significantly  different 
from  those  of  the  TM  and  TE  modes.  A  comparison  of  the  at¬ 
tenuation  rate  and  excitation  factor  modal  parameters  of  Mode 
X  with  those  of  Modes  1  through  4  has  revealed  that  the  con¬ 
tribution  of  the  Mode  X  signal  to  the  multimode  signal  is  sig¬ 
nificant  (and  must  be  included  in  the  multimode  signal  calcu¬ 
lations)  if  the  path  segment  is  in  night  and  has  a  westerly 
geomagnetic  bearing,  and  lies  in  the  equatorial  region.  This 
is  due  to  the  fact  that  the  "excitabi 1 ity/survivabi 1 ity "  of 
Mode  X  along  a  westerly-directed,  equatorial  path  segment  is 
comparable  to  that  of  Mode  1.  For  a  path  segment  that  is  not 
westerly-directed  and  equatorial,  attenuation  rate  of  Mode  X 
is  10  to  20  times  higher  than  that  of  Mode  1;  as  a  conse¬ 
quence,  Mode  X's  contribution  to  the  multimode  signal  along 
such  a  path  segment  would  be  insignificant  compared  to  that  of 
Mode  1,  and  hence  can  be  safely  ignored. 

33  SIGNIFICANT  SCATTERING  INTERFACE  TYPES 

This  section  presents  a  list  of  the  path-segment  in¬ 
terface  types  where  scattering  effects  must  be  computed  using 
a  mode  conversion  model.  For  these  interfaces,  the  eigen- 
value(s)  of  one  or  more  of  the  component  modes  of  the  multi- 
mode  signal  change  rapidly  across  the  interface,  and  thus 
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cause  significant  scattering  in  the  signal.  The  change  in 
eigenvalues  at  an  interface  is  related  to  the  change  in  path 
properties  at  the  interface.  The  eigenvalue/path-property 
sensitivity  plots  (Refs.  12  and  13)  can  be  examined  to  deter¬ 
mine  the  precise  segment  -  to  -  segment  change  in  eigenvalues 
along  any  portion  of  a  path.  This  examination  has  led  to 
identification  of  the  following  three  interface  types  where  a 
mode  conversion  model  must  be  used  to  compute  the  scattering 
effects : 

1)  Interface  at  the  day/night  terminator  cross¬ 
ing  along  a  path 

2)  Interface  in  the  equatorial  region  along 
westerly-directed  paths 

3)  Interface  with  a  discontinuity  in  the  path's 
ground  "conductivity  grade". 

The  ground  conductivity  grades  and  associated  ground  conduc¬ 
tivity  levels  are  given  in  Table  3.3-1.  In  this  table,  note 
that  conductivity  levels  5  through  10  are  grouped  together  and 
denoted  as  conductivity  grade  V  because  eigenvalues  of  the 
signal  component  modes  in  this  range  are  insensitive  to 
changes  in  the  conductivity  level  (Refs.  12  and  13). 

3.4  MSP  ALGORITHM  DEVELOPMENT 

This  section  presents  the  equations  used  by  the  MSP 
algorithm  for  computing  amplitude,  phase,  and  MIPD  of  an  Omega 
station  signal  for  a  given  location  and  time.  The  algorithm 
is  based  on  wavegu i de - mode  theory  with:  (1)  the  approximations 
and  assumptions  listed  in  Section  3.2,  (2)  the  significant 

.scattering  interface  types  listed  in  Section  3.3,  and  (3)  the 
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TABLE  3.3-1 

CONDUCTIVITY  GRADES  AND  ASSOCIATED  CONDUCTIVITY  LEVELS 


CONDUCTIVITY 

GRADE 

CONDUCTIVITY 

LEVELS 

CONDUCTIVITY 

VALUE 

(mho/m; 

I 

1 

10 '5 

I  I 

2 

3.2  x  10  "5 

I  I  I 

3 

10 

IV 

4 

3.2  x  10 '4 

V 

5  through  10 

10‘3  to  4 

computation  of  scattering  effects  outlined  earlier  in  Section 
3.1. 

Section  3.4.1  gives  a  list  of  the  symbols  used  in  the 
algorithm  equations.  Section  3.4.2  presents  the  equations  for 
computing  the  multimode  signal  along  a  path.  The  equations 
for  computing  the  amplitude,  phase,  and  MIPD  of  the  signal  are 
given  in  Section  3.4.3. 

3.4.1  List  of  Symbols 

This  section  presents  the  symbols  used  in  the  multi- 
mode  signal  prediction  algorithm  equations  contained  in  Sec- 
t ions  3.4.2  and  3.4.3. 
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S  vmbo  1 


Def init  ion 


E  (  s  ,  r  » 


Ek(s , r; 


e  r  (  s  > 


k ,  m  ,  n 


m 


'Ll’ 


m. 


Vertical  component  of  the  multimode  signal 
electric  field  at  path  segment  s  which  is 
at  a  distance  r  from  the  transmitter 

Vertical  component  of  the  kth  mode  signal 
electric  field  at  path  segment  s  which  is 
at  a  distance  r  from  the  transmitter 

Vertical  component  of  the  kth  mode  signal 
electric  field  at  the  point  just  before 
the  interface  between  path  segments  s  and 
(s+1) 

Mode  index  subscripts*  in  path  segments 
1,2,3,  respect  1  ve  ly 

Mode  index  subscripts  in  path  segments 
(N-l)  and  N,  respectively 


M  <  S  ) 


Number  of  modes  included  in  the  multimode 
signal  at  path  segment  s 


r  ( s  I 


Distance  f rom  transmitter  to  the  end  of 
path  segment  s 


Z 


Transmitting  antenna  constant 
(  =  i  2tt  Z  IG  h  Q  /  A. 3  2  j 

Free-space  characteristic  impedance 


Transmitting  antenna  current  amplitude 


'  . o  a vo  d  clutter  in  the  signal  equation,  the  mode  index  sub¬ 
scripts  k.  m ,  and  n  are  used  in  place  of  mi,  m2,  and  m3, 
respect : ve ly 
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Svmbol  icont 

ho 

A. 

a 

i 

s  in  1/2  (  r/a) 

B  v (  s  ) 


K 


3k(Sj 


a k ( s )  ,  yk  ( s ) 


i  E 

<t> 

4>. 


A*() 


Definition  ( cont . ) 

Effective  height  of  transmitting  antenna 
Signal  wavelength 
Earth's  radius 


V^T 


Geometrical  spreading  factor  describing  fo¬ 
cusing  of  signal  caused  by  spherical  nature 
of  the  earth 

Eigenvalue  (complex  quantity)  of  the  kth 
mode  at  path  segment  s 

{l-[H/2a] } ,  where  H  is  reference  height  at 
which  the  modified  refractive  index  of  the 
earth- ionosphere  waveguide  medium  is  unity 

[  7k(s)  +  i  ak(s)],  complex  value  of  the 
propagation  constant  of  the  kth  mode  at  path 
segment  s  at  the  earth's  surface,  as  de¬ 
fined  in  IPP;  it  is  a  function  of  the 
eigenvalue  of  the  kth  mode  at  path  segment  s 
and  the  path-segment  properties 

Attenuation  rate  and  phase  variation,  re¬ 
spectively,  of  the  kth  mode  at  path  segment 
s  on  the  earth's  surface,  as  defined  in  IPP 

Amplitude  of  the  complex  signal  field  E 

Phase  of  the  complex  signal  field  E 

Phase  of  the  complex  signal  field  Ei  of 
Mode  1  signal 

Modal  mt-irference-  induced  phase  deviation 
(MI  PD i  in  the  signal 
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Svmbo 1  ( com .  )  Definition  ( cont .  ) 

Smk(p)  Scattering  coefficient  at  interface  between 

segment  p  and  segment  (p  +  l)  for  the  kth  inci¬ 
dent  mode  and  the  mth  transmitted  mode;  it  is 
the  ratio  of  the  mth transmitted  mode  signal 
field  and  the  kth  incident  mode  signal  field 
at  the  interface;  it  is  computed  using 
either  FASTMC's  mode  conversion  model  or 
IPP's  WKB - approximat ion  depending  upon  the 
interface  type  (see  Section  3.3);  under  WKB - 
approximation:  Smk(p)  =Sm  k[  Xk  (p  +  l )  /Xk  (p)  l1/2, 
where  Sm  k  is  unity  if  m=k,  and  0.0  otherwise 

Xk(s)  Excitation  factor*  (complex  value)  of  the 

vertical  component  of  electric  field  of  the 
kth  mode  at  path  segment  s;  it  is  a  measure 
of  relative  efficiency  with  which  the  kth 
mode  is  excited  (or  received)  at  path  seg¬ 
ment  s  and  it  is  normalized  to  unity  for  a 
segment  with  infinite  ground  conductivity 
along  a  flat  earth 


'Note  that  the  excitation  factor  definition  in  the  MSP  algo¬ 
rithm  is  slightly  different  from  the  definition  used  in  IPP 
(Refs.  1  and  11).  The  difference  is  in  the  factor:  sin 
[9k(sj],  included  in  the  denominator  of  the  excitation  factor 
equation  (see  Refs.  1  and  14).  IPP  uses  a  9k(s)-value  de¬ 
termined  at  the  reference  height,  H,  of  70  km  (this  is  the 
height  where  the  modified  refractive  index  of  the  path  medium 
is  unity)  whereas  the  MSP  algorithm  uses  a  9k(s) -value  com¬ 
puted  at  the  earth's  surface.  Since  the  signal  field  of  in¬ 
terest  is  the  field  at  the  earth's  surface,  the  IPP-provided 
excitator  factor,  X£pp  (s),  is  inappropriate  for  use  in  the 
algorithm,  and  the  correct  value  is  that  value  used  in  the 
algorithm,  Xk(s;,  which  is  related  to  the  IPP-provided  value 
by  the  relationship:  Xk(s)  =  [X*?p  ( s )  ]  K3/2,  where  K  = 

{ 1- [  (H  =  70  j / ( 2a) ]  )  . 
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Svmbo I  fcont . ) 


Def init ion  (cont . ) 


Gk(s)  Height-gain  function*  (complex  value)  of  the 

vertical  component  of  the  electric  field  of 
the  kth  mode  at  path  segment  s,  as  defined  in 
IPP  (Refs.  1  and  11);  it  gives  the  variation 
of  the  electric  field  with  altitude  above 
the  earth's  surface  and  it  is  normalized  to 
unity  at  the  earth's  surface 


Note  that  both  the  height-gain  function,  Gk(s)  and  excita¬ 
tion  factor,  Xk(s)  are  defined  differently  in  IPP  (Refs.  1 
and  11)  and  FASTMC  (Ref.  9).  The  IPP-defined  and  FASTMC- 
defined  values  are  related  as: 


ipp 

Gv  (S) 


=  -(£ASTMC(s)  /sin  (  (6 k ( s )  ] 


X 


ipp 

k 

IF 
k 


(s)  =  -sin  [ 9  k ( s ) ]  X 


FASTMC, 


'(S) 


c.IPP  (S) 


xlPP  (S)  =  cfAS™c(s)  XPAS™C(S) 


Since  the  height-gain  function  and  excitation  factor  quan¬ 
tities  appear  as  a  product  in  the  MSP  equation  (see  Section 
3.4.3),  either  mode  1  - prov i ded  values  for  the  product  can  be 
correctly  used  in  the  signal  equation. 
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3.4.2  Algorithm  Equations 

Signal  propagation  equations  are  developed  below  for  a 
N-segment  long  path  shown  in  Fig.  3.4-1.  Each  path  segment  is 
characterized  by  constant  (homogeneous)  path  properties  associ¬ 
ated  with  the  mid-point  of  the  segment.  The  multimode  signal 
at  segment  s  is  assumed  to  be  composed  of  M(s)  modes.  For  the 
sake  of  clarity  in  presentation,  the  MSP  algorithm  equations 
for  the  path  shown  in  Fig.  3.4-1  are  presented  below:  first, 
along  path  segment  1;  second,  along  path  segment  2;  third, 
along  path  segment  3;  and  fourth,  along  path  segment  N. 


G-0VI2" 
4-  •  5  ' 


Figure  3.4-1  \-Segment  Path  Displaying  Mode  Conversion/ 

Scattering  at  Path-Segment  Interfaces  Along  a 
N-Segment  Long  Path 
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Note  that  the  excitation  factor,  Xk(s),  in  Equation 

3.4-1  and  elsewhere  in  this  report  is  the  IPP-provided  value 

5/2 

multiplied  by  K  (see  footnotes  to  the  excitation  f  CL  C  o  O  i  a  ltd 

height-gain  function  definitions  presented  in  Section  3.4.1). 


Signal  Field  Along  Path  Segment  1 

The  multimode  signal  field  along  path  segment  1 
(which  includes  the  transmitter  but  exclude*  the  interface  be¬ 
tween  path  segments  1  and  2),  at  a  distance  r  (o  <  r  <  r  (lj) 
from  the  transmitter  along  the  path  shown  in  Fig.  3.4-1,  is 
given  by: 


y.D 

E(k,r)  =  £  Ek(l,r) 
k  =  l 


sin 


-1/2 


(r/ah 


Freq.  h 
Trans . 
Power  - 
Dependent 
Constant 


Geometric- 

Spreading 

Factor 


M  ( 1 ) 

I 


k=l 


Trans  .  Rec  . 


Exc i tat i on 
Factors 


Trans  .  Rec  . 


Height-Gain 
Funct ions 


exp  |  i. (3 k ( 1 )  [r]  | 


Propagation 
Factor  of  kth  mode 
in  Segment  1  over 
distance  r 


1/2 


(3.4-1) 


Signal  Field  Along  Path  Segment  2 

The  multimode  signal  field  along  path  segment  2  (in¬ 
cluding  the  interface  between  path  segments  1  and  2  but  ex- 
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eluding  the  interface  between  path  segments  2  and  3),  at  a 
distance  r  (r  (1)  <  r  <  r  (2)),  from  the  transmitter  along  the 
path  shown  in  Fig.  3.4-1,  is  given  by  extending  Eq.  3.4-1  as 
described  below: 


M(  2) 

E(2,r)  =  Y  E  (  2  ,  r ) 
m 

m=l 

where 


E 


m 


(2  ,  r) 


M(  1 ) 

I 


k=l 


E 

kth  mode 
field  just 
before 
interface 
between 
Segments 
1  &  2 


Scatt.  coeff. 
at  interface 
between  segments 
1  &  2  describing 
conversion  of 
kth incident 
mode  signal 
into  mth trans¬ 
mitted  mode 
signal 


(3 . 4 -2a) 

[  G.(2)/Gk(l)] 

'*■ . .  ™  ™  V 

Rat io  of  height  - 
gain  functions 
to  account  for 
change  in  G  from 
segment  l  to  seg¬ 
ment  2 


exp|i3m  (  2)  [  r-r  ( 1)  ]  J-  sin  1/2  [  r/aj/s  in 


-1/2 


[  r ( 1 ) /a] 


} 


Propagation  Factor 
of  mthmode  over 
(r-r(l)]  distance 
in  Segment  2 


Ratio  of  Spreading  Factors  to 
account  for  change  in  the  station 
distance  from  segment  l  to  seg¬ 
ment  2 


(3 . 4-2b) 


There  f ore  , 

E  i  2  ,  r  j  =  CQ  s  irr1/2  (  r/a) 


X 


m=  1 


X 


m ,  k 


(1) 


GJ2; 


mu)  r 

i  [ 

k=l  L 


X,  (1)  G,  (1)  exp  4  ip,  (1)  [r(l) 


{ 


exp  ^  i 3 m ( 2 )  [ r-r ( 1) 


{ 
}] 


>} 


(3.4-2C) 


3.4-2 


Note  that  the  scattering  coefficient  S  ,  (1)  in 

t  h  m  ,  R 

describes  the  conversion  of  the  k  incident  mode 


Eqs  . 
s  igna  1 
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field  into  the  transmitted  mode  signal  field  due  to  mode 

conversion  at  the  interface  between  path  segments  1  and  2. 

The  scattering  (or  mode  conversion)  is  due  to  a  step  change  in 
the  path  properties  at  the  interface.  Also  note  that  for  the 
interface  with  no  discontinuity  in  the  mode-eigenvalue, 

S  ,  (l)  =  8  ,  ,  where  8  ,  is  Kronecker  delta  function  which 

is  unity  for  m=k  and  zero  for  m^k.  Depending  upon  the  extent 

of  mode - e igenva lue  discontinuity  at  an  interface,  the  scatter¬ 
ing  coefficients  are  computed  using  the  FASTMC's  mode  conver¬ 
sion  model  (Ref.  9)  or  IPP's  WKB - approx imat ion  approach  (Ref. 
11),  as  discussed  in  Sections  3.3  and  3.4.1. 

Signal  Field  in  Path  Segment  3 

The  multimode  signal  field  along  path  segment  3  (in¬ 
cluding  the  interface  between  path  segments  2  and  3  but  ex¬ 
cluding  the  interface  between  path  segments  3  and  4),  at  a 

distance  r  (r  (2)  <  r  <  (3))  from  the  transmitter  along  the 
path  shown  in  Fig.  3.4-1,  is  given  by  extending  Eqs .  3.4-2  as 
described  below.- 


M  ( 3  ) 

E  (  3  ,  r  )  =  Y  E  (  3  ,  r ) 
1  n 
n=l 


(  3 . 4  -  3  a  ) 


where 


En(3,r) 


Mr  2> 
V 

z, 

m=l 


ET2>  [<=„  <3>/Gm(3,  ] 


exp-fi(3ri(3}[r-r.2))J.-|sin  1/2[r/a|/sin  1/2[r(2)/a;| 


(3 .4 -3b) 
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Therefore , 


E(3,r)  =  CQ  sin'1/2(r/a) 

M(3) 

X 

M  ( 2 ) 

I 

M(l) 

X 

xk(D  G  k  ( 1 ) 

n=l 

m=l 

k  =  l 

X 


X 


exp|i(3k(l)  { r ( l )  ]| 

Sn.m(2>  Gn<3>  e*P 


sm  k(i)  exp|i3m(2)  [r(2)-r(l)  ]} 
|ipn(3)  [r-r(2)  ]  J-l  (3.4-3' 


Signal  Field  in  Path  Segment  N 

The  multimode  signal  field  along  path  segment  N  (in¬ 
cluding  the  interface  between  path  segments  (N-i)  and  N,  but 
excluding  interface  between  path  segments  N  and  (N+l)),  at  a 


distance  r 

(r(N-l) 

<  r 

path  shown 

in  Fig. 

3.4- 

described 

below 

M(N) 

E(N,  r)  =CQ  sin'1/2(r/a) 

mN  =  1 

x  • 

exp 

X  - 

exp 

x  \ 

[  S  (  N—  1 ) 

G(N) 

m,, 

|xk(l)  Gk(l)  exp[i3k(l)  (r(l)  f]} 


t 

exp  (N)  [r-r(N-l)  ]  jJ- 


(3.4-4) 
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In  Eq.  3.4-4,  nr,  and  m..  .  denote  the  mode  indices  of 

N  N-l 

the  component  modes  of  the  signal  along  path  segments  N  and 
(N-l),  respectively;  S  (N-l)  is  the  scattering  coeffi- 

N  N-l 

cient  at  the  interface  between  path  segments  (N-l)  and  N  for 

incident  and  transmitted  modes  at  the  interface. 

N-l  N 


3.4.3  Signa 1 


MI  PD 


iiat  ion; 


This  section  develops  equations  for  computing  the  am¬ 
plitude,  phase,  and  MIPD  (modal  interference  -  induced  phase  de¬ 
viation)  of  the  multimode  signal  from  the  signal  field  equa¬ 
tions  given  in  Section  3.4.2. 


Denoting  the  in-phase  and  quadrature  components  of  a 
signal  field  by  the  superscripts  i  and  q,  respectively,  the 
amplitude,  |E|,  and  phase,  <j>,  of  the  complex  signal  field,  E, 
are  given  by: 

|E|  (Amplitude  of  E)  *  (E2  +  E2)172  (3.4-5a) 

4>  (Phase  of  E)  =  tan  ~1  ( Eq /E 1  )  (3.4-5b) 

Thus,  the  phase  of  the  Mode  1  component  signal  field 
E1  of  the  multimode  signal  E  is  given  by: 

$  =  tan"1 (E^/Ej )  (3.4-6) 

The  in-phase  and  quadrature  components  of  E  and  E^ 
are  computed  by  summing  over  the  component  modes: 


E  s  in  6 
p  p 


(3.4-7) 
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E*  =  Ex  COS  4>x 

E?  =  E,  sin  <j> 

-L  1  1 


(3.4-8) 


where  P  is  the  number  of  component  modes  of  the  multimode 


s igna. 


Thus  MI  PD ,  A(f>,  is  given  by 


A<t>  =  4>  -  4>1 


(3.4-9) 


Figure  3.4-2  shows  MIPD  for  the  case  where  the  multi- 
mode  signal  is  composed  of  two  modes  (Modes  1  and  2).  For  the 
case  when  the  multimode  signal  is  composed  of  more  than  two 
modes  (e.g.,  P  modes),  Fig.  3.4-2  can  be  used  to  determine  MIPD 
by  assuming  that  E2  is  the  phasor-sum  of  Modes  2  through  P. 


Figure  3.4-2  Multimode  Signal  E  as  a  Phasor-Sum  of  Ex  and 

Ex  (Modes  1  and  2) 
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4  . 


SUMMARY 


Omega  navigation  is  based  on  the  assumption  that  an 
Omega  signal,  theoretically  comprised  of  an  infinite  number  of 
modes  (in  the  full-wave  waveguide -mode  theory  representation), 
is  adequately  approximated  by  the  signal's  Mode  1  component. 
This  is  valid  at  a  majority  of  geographic  locations  and  times. 
The  lack  of  Mode  1  dominance  in  a  signal  is  known  as  "modal 
interference"  and  the  resulting  deviation  of  the  multimode 
signal  phase  from  the  signal's  Mode  1  component  phase  is 
called  modal  interf erence - induced  phase  deviation  (MIPD) . 
Successful  use  of  an  Omega  navigation  signal  requires  that  the 
signal  be  strong  enough  such  that  its  SNR  exceeds  the  receiver 
detection  threshold  level  and  the  signal's  MIPD  be  below  a 
certain  threshold.  The  ONSC- recommended  value  of  that  thresh¬ 
old  is  20  centicycles. 

The  Analytic  Sciences  Corporation  is  developing  a 
worldwide  database  for  the  Omega  Navigation  System  Center, 
consisting  of  the  model -predicted  values  of  the  amplitude  and 
MIPD  of  10.2  and  13.6  kHz  Omega  station  signals  for  96  repre¬ 
sentative  times  of  the  year.  The  database  can  be  developed 
from  the  Naval  Ocean  System  Center's  FASTMC  computer  model 
(Ref.  9)  but  this  approach  is  computationally  very  expensive. 
The  objective  of  the  effort  reported  herein  is  to  develop  a 
cost-effective  Multimode  Signal  Prediction  (MSP)  algorithm 
(based  on  waveguide-mode  theory)  for  computing  the  amplitude 
and  MIPD  of  an  Omega  station  signal  for  a  given  location  and 
time. 
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Since  an  Omega  signal  path  is  i nhomogeneous ,  the  al¬ 
gorithm  approximates  the  path  as  a  concatenation  of  a  large 
number  of  homogeneous  path  segments  each  with  constant  elec¬ 
tromagnetic  properties  associated  with  the  mid-point  of  the 
segment.  The  signal  along  a  day/night  path  segment  is  ap¬ 
proximated  in  the  algorithm  by  the  sum  of  the: 


•  Three  lowest  phase- velocity  modes  (Modes  1 
through  3)  of  the  daytime  earth- ionosphere 
waveguide  modeling  the  day  path  segment 

•  Five  lowest  phase-velocity  modes  (Modes  1 
through  4,  and  Mode  X*)  of  the  nighttime 
earth- ionosphere  waveguide  modeling  the 
night  path  segment. 


The  MSP  algorithm  consists  of: 


(1)  Determining  modal  parameters  of  the  selected 
component  modes  of  the  multimode  signal  for 
each  of  the  homogeneous  path  segments  of  the 
signal  path;  these  are  conveniently  obtained 
using  the  TASC-developed  look-up  table  of  the 
multimode  signal  modal  parameters 

(2)  Computing  scattering  effects  caused  by  mode 
conversion  at  each  of  the  interfaces  (between 
adjoining  path  segments)  along  the  signal  path 
using:  (a)  the  FASTMC  mode  conversion  model 
(Ref.  9)  if  the  eigenvalue  of  any  one  of  the 
component  modes  of  the  signal  varies  suffi¬ 
ciently  rapidly  across  the  interface,  other¬ 
wise  (b)  using  the  computat ional ly - e f f ic ient 
WKB - approx imat ion  used  by  IPP  (Ref.  10) 

■*  Mode  X  is  a  new  mode  recently  identified  by  TASC  (Refs.  12 
and  13).  The  propagation  characteristics  of  this  mode  are 
comparable  to  that  of  Mode  1  along  westerly-directed  path- 
segments  (with  geomagnetic  bearings  between  190°  and  350°, 
measured  clockwise  from  the  geomagnetic  north)  in  the  equa¬ 
torial  region  (between  10°S  and  10°N  geomagnetic  latitudes) 
Therefore,  Mode  X  must  be  included  in  the  multimode  signal 
calculations  for  these  path  segments;  elsewhere  it  can  be 
ignored . 
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(3)  Modifying  the  individual  path-segment  multi- 
mode  signal  propagation  characteristics  ob¬ 
tained  in  step  (1)  by  incorporat ing  the 
scattering  effects  computed  in  step  (2)  at 
the  segment  interfaces. 

The  resulting  MSP  algorithm  is  computationally  much  more  effi¬ 
cient  than  the  FASTMC  computer  model  (Ref.  9).  This  is  because 
the  MSP  algorithm  computes  scattering  effects  at  a  path-segment 
interface  using  the  mode  conversion  model  only  if,  based  on  a 
priori  mode - e igenvalue  variation  information  across  the  inter¬ 
face,  the  effects  are  expected  to  be  large;  otherwise  it  uses 
the  WKB-approximation  of  IPP  (Ref.  11).  In  nearly  all  cases, 
the  algorithm  has  been  found  to  provide  signal  predictions 
within  one  dB  of  the  FASTMC  predictions  at  about  one  tenth  of 
the  FASTMC 's  computer  cost. 
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